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Abstract 

After reviewing the motivation for high-statistics charm studies, 
we describe a fixed-target experiment capable of reconstructing > 10^ 
charm decays. High-rate sihcon and scintillating-fiber tracking sys- 
tems allow operation at a 5 MHz interaction rate, using an 800 GeV 
primary proton beam incident on a thin target. The trigger includes 
transverse energy and/or lepton requirements and also requires evi- 
dence of decay vertices. Simulations show that adequate rejection can 
be obtained. 



1 Charm Physics at Fermilab: Brief History 
and Outlook 

The program of fixed-target charm studies at Fermilab which began in earnest 
about 1984 with E691 has been by any measure extraordinarily successful. 
The history of Fermilab charm experiments is illustrated in Fig. |l| which 

*Presented at the LISHEP95 Workshop, Rio de Janeiro, Brazil, February 20-22, 1995. 
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shows roughly exponential growth in sensitivity since the late 1970s. While 
the physics reach of such experiments depends both on the number of signal 
events reconstructed and on the amount of background under the peaks, the 
former figure can still serve as a starting point for discussion. This number 
is expected to reach ~ 10^ events during the next few years with the runs of 
Fermilab E781 and E831 (and the advent of CLEO III). 
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Figure 1: Yield of reconstructed charm vs. year of run for completed and 
approved high-statistics Fermilab fixed-target charm experiments; symbols 
indicate type of beam employed. 

We are naturally led to the question whether further substantial advances 
in sensitivity can be foreseen in the years beyond. Our conclusion (detailed 
below) is that an optimized fixed-target apparatus can reach ~ 10^ recon- 
structed decays using technology which is now or will soon be available. At 
these levels of sensitivity charm studies will provide increasingly incisive tests 
of the Standard Model and will also begin to probe significantly beyond it [@]. 



2 



1.1 Searches for new physics 

Given the large number of free parameters in the Standard Model, we must 
suspect that some new physics lies beyond it. This new physics could have de- 
tectable effects in the charm sector. For example, Pakvasa 0] has emphasized 
that rare and forbidden processes such as — mixing, charm-changing 
neutral currents, and even lepton-family-violating currents must occur at 
some level if we are ever to understand the pattern of fermion masses and 
mixings. While direct CP violation in singly- Cabibbo-suppressed (SCSD) 
charm decays is expected ||^ in the Standard Model at the ^ 10~^ level (pos- 
sibly observable at the level of sensitivity discussed here), CP violation in 
Cabibbo-favored or doubly-Cabibbo-suppressed modes would be a clear sig- 
nature of new physics. 

We can identify a growing trend towards precision studies of the strange 
quark (e.g. E871 and the KTeV and KAMI programs at Fermilab and E787 
and E865 at the Brookhaven AGS) and the beauty quark (CDF and DO 
upgrades, HERA-i?, and the SLAG and KEK B Factories). The major goal 
of these efforts is the detailed study of suppressed Standard-Model effects 
{CP violation and rare decays), where small deviations may reveal aspects 
of the ultrahigh-energy theory. Various proposed extensions of the Standard 
Model (e.g. multiple-Higgs-doublet models with [|, |[ or without tree- 
level flavor-changing couplings, left-right-symmetric models f^, and models 
with family symmetry or supersymmetry |T^) predict effects also in the 



charm sector which should be detectable at the level of sensitivity discussed 
here. Gharm experiments are uniquely sensitive to couplings of these new 
currents to "up-like" quarks, which may differ from those to "down-like" 
quarks |]n|. Furthermore, the Standard-Model "backgrounds" in processes 
such as particle-antiparticle mixing, flavor-changing neutral currents, and CP 
violation are far smaller in charm than in strangeness and beauty, resulting 
in a larger and different discovery window. Precision studies of charm are 
thus complementary to those of strangeness and beauty. 

jjo _ 1)0 ]2iixing may be one of the more promising places to look for 
low-energy manifestations of physics beyond the Standard Model. Standard- 
Model contributions to |Am£)| are estimated 0] to give rmix < 10~^. Many 
nonstandard models predict a much larger effect. An interesting example 
is the multiple-Higgs model of Hall and Weinberg 0, in which |Am£)| can 
be as large as 10~^eV, approaching the current experimental limit. In this 
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model CP violation arises from the Higgs sector, and CP violation in the 
beauty sector is expected to be small. 

Charm-changing neutral currents are forbidden at tree-level in the Stan- 
dard Model and thus are predicted [|12| to have branching ratios of order 
10~^^ to lO"*^. A variety of non-Standard models [0 predict effects sub- 
stantially larger than this. Experimental sensitivities are now in the range 
~ 10~^ and are expected to reach ~ 10~^ in E831. 

In summary (based on discussions and presentations at the CHARM2000 
Workshop [|14|), we anticipate the following order-of-magnitude estimates for 
sensitivity to new physics in a "10^-charm" experiment: 

• mixing: rmix ~ 10^^ 

• charm-changing neutral currents and lepton-family- violating currents: 
B ~ 10"^ 

• charm CP violation: Aqp ~ 10~^ in SCSD modes 



1.2 Tests of the Standard Model 

Incisive tests of the CKM Model are envisioned from precision beauty studies 
at CLEO, the Tevatron Collider, HERA-i?, and the B factories. Here too, 
high-sensitivity charm-quark studies play an important role complementary 
to that of b experiments: 



1. Several authors |]I5|, |I6[ have noted that suppressed Standard-Model 
processes (including such nonperturbative effects as final-state inter- 
actions) can complicate attempts to extract CKM parameters from 
beauty measurements. Charm decay provides an ideal laboratory in 
which to study these effects, which are both larger (due to the smaller 
charm-quark mass) and more cleanly studied (due to the small size of 
Standard-Model CP-violating and penguin-mediated effects) in charm 
than in beauty. In addition, the interplay of perturbative and nonper- 
turbative Standard- Mo del physics at work here is of interest in its own 
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2. The extraction of CKM parameters from beauty measurements in many 
cases requires knowledge of form factors and decay constants. These 
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have been predicted using various models, including lattice gauge calcu- 
lations and calculations employing the heavy-quark symmetry of QCD. 
The predictions of these models need to be tested in the charm sector 
to verify their applicability in the beauty sector. 

3. Charm is the lightest quark for which perturbative QCD calculations 
may be meaningful, thus its production in large numbers in hadronic 
interactions can provide detailed tests of QCD. Understanding the pro- 
duction mechanism in terms of perturbative-QCD processes could lead 
to a determination of the gluon distribution of the proton at x 0.025, 
complementary to the small-x measurements possible at HERA and 
the hadron-coUider experiments. The interplay of perturbative and 
non-perturbative effects here presents the opportunity to study longer- 
range aspects of QCD, such as the potential existence (and magnitude) 
of an intrinsic charm component in the nucleon, and final-state inter- 
actions in propagation of charm quarks through nuclear matter. 



2 A Next-Generation Charm Spectrometer 

We next discuss a hypothetical experiment which potentially can achieve 
the 10^-reconstructed-charm sensitivity mentioned above. As we will see, 
the most demanding requirement is on the trigger. In particular, an on- 
line secondary-vertex trigger is needed if adequate trigger rejection is to be 
achieved without sacrificing sensitivity in hadronic decay modes. 



2.1 Beam and target 

To achieve the sensitivity discussed here in a fixed-target run of ~ 10^ 
beam spills requires a primary proton beam (see Table 1 ll8|). Given 
a{pN DX) + a{pN -> D X) ^ 40 /ib/nucleon at 800 GeV and 



a {pN — > X) oc A^-^ [^, and assuming that the cross section to produce 



Ds and charmed baryons is ~ 15% that of D mesons, the charmed-particle 
production rate is 7 x 10~^/interaction if a high-A target (e.g. Au) is used. 

A target which is short compared to typical charm decay lengths is crucial 
for optimizing background suppression, both off-line and at trigger level. 
While multiple thin targets could be employed (as in E791 and E831), a 
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Table 1: Charm yields* for various possible beams using an open-geometry 
general-purpose spectrometer. 



Beam 


0"cc 


Max. beam 


Charm events 


Reconstructable 




(cm^) 


/ spill 


/ spill 


charm/spill^ 


p (800 GeV) 


2 X 10-29 


10^3 


7x10^^ 


^ 1.5 X 10^^ 


vr- (500 GeV) 


2 X 10-29 


108 § 


3 X 10^ 


^ 60 


7 (^200 GeV) 


1 X 10-3° 


108 § 


2 X 10^ 


^ 4 



* neglecting A-dependence enhancement. 

^ based on E (branching ratio x efficiency) w IQ-^ for reconstructable modes. 

^ assuming a maximum of 0.1 interaction/bucket. 

S limited by proton economics. 



single target facilitates fast vertex triggering. A 1mm W, Pt, or Au target 
is one possibility, representing Ril% of an interaction length and on average 
14% of a radiation length for outgoing secondaries. A low-Z material such as 
i^C-diamond may be favored to minimize scattering of low-momentum pions 
from D* decay |TB[; then a 2 mm target is suitable, representing ^1% of an 
interaction length and ~1% of a radiation length and producing charm at 
the rate 3 x lO-^/interaction. Given the typical Lorentz boost 7 ~ 20, a 
1-2 mm target is short enough that a substantial fraction even of charmed 
baryons will decay outside it. 

We consider a benchmark 5 MHz interaction rate, which then requires 
0.5-1 GHz of primary proton beam, an intensity easily attainable. As shown 
below, this yields ^ 10^ reconstructed charm per few x 10^ s of beam (^ 10^ 
spills X 20s/spill). 



2.2 Rate capability 

A significant design challenge is posed by radiation damage to the silicon 
detectors. To configure detectors which can survive at the desired sensitiv- 
ity, we choose suitable maximum and (in one view) minimum angles for the 
instrumented aperture, arranging the detectors along the beam axis with a 
small gap through which pass the uninteracted beam and secondaries below 
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the minimum angle (Figs. 0, ^).[] Thus the rate is spread approximately 
equally over several detector planes, with large-angle secondaries measured 
close to the target and small-angle secondaries farther downstream. Along 
the beam axis the spacing of detectors increases geometrically, making the 
lever arm for vertex reconstruction independent of production angle. Since 
small-angle secondaries tend to have high momentum, the multiple-scattering 
contribution to vertex resolution is also approximately independent of pro- 
duction angle. We have chosen an instrumented angular range < 200 mr, 
4 < 6y < 175 mr, corresponding to the center-of-mass rapidity range \y\ ^ 1.9 
and containing over 90% of produced secondaries. 
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Figure 2: Spectrometer layout (bend view) 



To maximize the rate capability of the spectrometer, tracking is per- 
formed entirely with silicon and scintillating-fiber planes. The rate per unit 
area (and hence the radiation fluence) in a detector element can easily be 
estimated based on the uniform-pseudorapidity approximation. Fig. ^ shows 
the rate calculation for an annular area dA located a transverse distance r 
from the beam and of thickness dr. Since the operational limit of present- 
day silicon detectors is 10^^ particles/cm^, the charged multiplicity per unit 
pseudorapidity in 800 GeV proton-nucleus collisions is n ~ 4 for high-A tar- 
gets (less for C), and calling the total number of interactions in the run 



^An alternative approach with no gap may also be workable if the beam is spread over 
sufficient area to satisfy rate and radiation-damage limits, however the approach described 
here probably allows smaller silicon detectors and is "cleaner" in that the beam passes 
through a minimum of material. 
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Figure 3: Detail of vertex region (showing optional optical impact-parameter 
trigger). 



^int, we can derive the "minimum survivable" inner detector radius 
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Since a typical run will yield fewer than (5 x 10^ interactions/ s) x (4 x 10^ s) = 
2 X 10^^ interactions, we obtain conservatively r^^^^ = 3.5 mm, which we set as 
the half-gap between the two detector arms. This ensures that the detectors 
will survive for the entire run (or at most will need to be replaced onceQ). 
To cover the desired angular range, we configure 14 double-sided silicon strip 
detector^ above and 14 below the beam as shown in Fig. ^, such that at all 
angles of interest there are at least six measurements per track (and more at 
small angles where the occupancy is highest). 

A large-scale scintillating-fiber tracker has recently been successfully 
tested P2[ for the DO upgrade, and we base our fiber-tracker design on that 



demonstrated technology. Green-scintillating 3HF/PTP fibers are deployed 
in staggered doublets in three views. They are read out using cryogenic solid- 
state "visible-light photon counters" (VLPCs) [Q. Since the fibers are more 
radiation-hard than silicon detectors, and (due to occupancy; see below) the 

^In E789 we operated silicon detectors at fluence up to ~ 5 x 10^^ cni^^ with negligible 
efficiency loss. 

■^We assume silicon strip detectors for definiteness, but silicon pixel detectors would be 
better if available with sufficient readout speed and radiation hardness; because of their 
radiation hardness, diamond detectors |24], if available, should also be considered. 
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Figure 4: Calculation of rate per unit area in an annulus. 



beam gap between fiber planes is larger than that in the silicon, radiation 
damage of the fibers is not anticipated to be a problem. 

We assume 1-bucket (<19ns) recovery times for all detectors, so that 
there is no pile-up due to out-of-time interactions. Designs capable of this 
performance have been presented [21, 25] for all detectors except the TRD.0 

Detector-element occupancies also follow from the derivation of Fig. ^ 
For an element of height dy located a transverse distance y from the beam and 
covering — Xmax < x < x^ax, the occupancy per event (neglecting magnetic 
bending) is 

^ dy -^max 

arctan . 

Ti y y 

For 800 /im fiber diameter, this implies ~16% occupancy at y = 1 cm, ^8% 
at 2 cm, and ~4% at 4 cm. A full trackfinding simulation will be required 
to assess the maximum acceptable occupancy, but this suggests ~1 cm as 
the minimum acceptable half-gap in the scintillating-fiber planes. The fibers 
near the gap could be split at x = and read out at both ends, halving their 
occupancies. Since shorter fibers have less attenuation, a smaller diameter 
could be used near the gap, reducing occupancy still further. 



2.3 Spectrometer performance 

We have carried out a simple Monte Carlo simulation of the spectrometer 
sketched above. Assuming a 1.2-m-long analyzing magnet with pole pieces 

''it may be that a TRD for electron identification is not cost-effective and a hadron-blind 
detector |p6| or preshower detector should be used instead. 
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tapered to give 0.5 GeV pt kick, we obtain (56 ± 1)% geometrical acceptance 
for D° K-TI+ decays and (44 ± 1)% for D*+ D^Tr+ K~Tr+n+, com- 
parable to those of existing open-geometry spectrometers despite the beam 
gap. With silicon detectors of 25 /im pitch read out digitally (i.e. no pulse- 
height information) and 800 /xm scintillating-fiber pitch, and assuming ±10° 
stereo, Gaussian fits to the reconstructed distributions give rms resolutions of 
6 MeV in mass (a factor ^ 2 better than that of existing spectrometers) and 
11 /xm (bend-view) and 21 /xm (nonbend-view) in impact parameter, giving 
40 fs decay proper-time resolution, comparable to that of existing spectrom- 
eters. Since the mass resolution is dominated by scattering, minimization of 
material is crucial, for example use of helium bags and avoidance of thresh- 
old Cherenkov counters employing heavy gas mixtures. (The performance 
parameters just given are a snapshot of work in progress and probably can 
be improved with further optimization.) 



2.4 Trigger 

While the most successful previous charm hadroproduction experiments 
(E769 and E791) used very loose triggers and recorded most inelastic interac- 
tions, this approach is unlikely to extrapolate successfully by three orders of 
magnitude! (Consider that E791 recorded 2 x 10^° events - tens of terabytes 
of data - on 20,000 8 mm tapes.) Thus our sensitivity goal requires a highly 
selective trigger. However, we wish to trigger on charm-event characteristics 
which bias the physics as little as possible. Lepton triggers, used successfully 
by E653, while capable of great selectivity (~ 10^ rejection for minimum- 
bias events), have only ~ 10% charm efficiency. The Et triggers used by 
E769 and E791, while highly efficient for charm, have poor selectivity (^ 10 
minimum-bias rejection). We therefore assume a first-level trigger requiring 
calorimetric Et OR'ed with high-p^-lepton and lepton-pair triggers. At sec- 
ond level, secondary- vertex requirements can be imposed on the i?j-triggered 
events to achieve a rate (~ 100 kHz) which is practical to record. 

Analyses of the efficacy of an Et trigger carried out using E791 data 



28| and the PYTHIA Monte Carlo [^] agree on minimum-bias rejection vs. 



charm efficiency (though due to nuclear effects not simulated in PYTHIA, 
they differ as to the Ef threshold corresponding to a given rejection). Fig. ^ 
shows the efficiencies for charm and minimum-bias events as a function of 
the PYTHIA Et threshold. A considerable degradation results if there is 
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significant probability for two interactions to pile up in the calorimeter. Given 
the 53 MHz rf structure of the Tevatron beam and the typical ~ 50% effective 
spill duty factor, at the benchmark 5 MHz mean interaction rate there is a 
~ 20% probability for a second simultaneous interaction. Thus at a 5 GeV 
PYTHIA Et threshold (corresponding to a :^10GeV actual threshold), the 
minimum-bias rejection factor is 5, i.e. pile-up degrades the rejection by a 
factor ^ 2, even for a calorimeter with one-bucket resolution. The charm 
efficiency at this threshold is about 50%, for a charm enrichment of ~ 2.5. 
(These are rough estimates based on a relatively crude calorimeter and 
an optimized calorimeter may provide better rejection.) Such an Et trigger 
yields a 1 MHz input rate to the next level. 
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Figure 5: Minimum-bias and D efficiencies vs. Et threshold; points are from 
PYTHIA simulation of Ref. ||29|, curves are fits of the form exp {—aE^), and 
arrows indicate the 5 GeV "PYTHIA Et" threshold discussed in text. 



While it may be technically feasible by the Year 2000 to record events 
at a IMHz rate, an additional factor ^10 in trigger rejection is desirable, 
and can be achieved by requiring evidence of secondary vertices. Existing 
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custom trackfinding trigger processors while perhaps capable of this 
rejection, typically fall short by ~ one order of magnitude in speed. At 
~ 1 MlPS-s/event, an on-line farm of commercial processors would need a 
capacity of ~ 10^ — 10^ MIPS, which may be prohibitive even in the Year 
2000. It is likely that by then a sufficiently fast custom trackfinding processor 
can be developed. This would require fast buffering (~ 100 ns) and readout 
(~ 1 /is) of event information in order not to impose excessive deadtime. 
Trackfinding secondary-vertex triggers benefit from the use of focused beam 
and a single thin target, which allow simplification of the algorithm since 
the primary vertex location is known a priori. Since low-pt tracks have 
poor vertex resolution |^ , a trigger which discriminates pt is more effective 
than one which is purely topological; such discrimination may be simply 
accomplished by placing the vertex detectors in a weak magnetic field and 
looking for straight tracks. Christian has suggested a simple trigger- 
processor algorithm based on this idea. 

As a first step in evaluating the efficacy of such a trigger, a Monte Carlo 
simulation was performed for the silicon-detector configuration shown in 
Fig. ^ using PYTHIA for the event generation. The target length was taken 
to be 1 mm, while its transverse dimensions were ignored, approximating a 
very narrow target (at least in one view). 25 /xm strip pitch was assumed. 
An algorithm was applied that searched for hits in two seed planes in the y-z 
view. When a pair of hits was found with more than some minimum impact 
parameter when extrapolated to the target, additional y hits, consistent with 
the line defined by the two seed planes, were sought. The search was repeated 
using different pairs of seed planes, and an event was accepted if at least one 
track was found. Fig. ^ shows the efficiency of the algorithm, both for charm 
and for minimum-bias events, as a function of the minimum impact param- 
eter. It can be seen that for an impact-parameter cut of 200 /im, a rejection 
factor of more than 20 can be obtained for the background, while retaining 
a 67% efficiency for charm events. Somewhat surprisingly, the addition of a 
weak magnetic field (up to 3 kG) in the target region was found not to be 
helpful: because the algorithm looks for straight tracks, a field strong enough 
to improve the rejection degrades the efficiency significantly. 

A key feature of Christian's algorithm is the elimination of hits which lie 
on straight lines pointing to the target prior to searching for tracks of finite 
impact parameter. This can reduce substantially the processing time, since 
as the number of hits per detector plane {n) increases, the time to eliminate 
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Figure 6: Efficiency of the vertex-trigger algorithm as a function of the min- 
imum impact parameter for charm events (dots) and minimum-bias events 
(triangles), using a PYTHIA-based Monte Carlo simulation. 



hits is linear in n, while the time to find tracks of finite impact parameter 
goes as Ti? (due to the loops over hits in two seed planes) The results 
using this preliminary hit-elimination pass are shown in Fig. |^. For the 
200 /im impact-parameter cut discussed above, the rejection and efficiency 
are hardly affected. This algorithm thus looks quite promising. 

As alternatives to iterative trackfinding at a 1 MHz event rate, three other 
approaches also appear worth pursuing. The ffist is a secondary-vertex trig- 
ger implemented using fast parallel logic, e.g. PALs, neural networks, or pre- 
downloaded fast RAMs, to look quickly for patterns in the silicon detectors 
corresponding to tracks originating downstream of the target. The others 
are fast secondary-vertex trigger devices originally proposed for beauty: the 
optical impact-parameter and Cherenkov multiplicity-jump I^S] triggers; 
while results from prototype tests so far suggest undesirably low charm ef- 
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Figure 7: Efficiency of the vertex-trigger algorithm as a function of the min- 
imum impact parameter for charm events (dots) and minimum-bias events 
(triangles), using a preliminary hit-elimination pass to speed up processing. 



ficiency, these might with further development provide sufficient resolution 
to trigger efficiently on charm. For example, Fig. |^ shows the efficiency for 
minimum-bias, charm, and beauty events projected for a version of the opti- 
cal trigger jS^ , indicating 40% charm efficiency for a factor 5 minimum-bias 
rejection. The resulting ^200 kHz event rate can be processed on-line or 
recorded using existing technology. 



3 Yield 

The charm yield is straightforwardly estimated. Assuming a Au target and a 
typical fixed-target run of 3 x 10^ live beam seconds, 10^^ charmed particles 
are produced. The reconstructed-event yields in representative modes are 
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Figure 8: Estimated optical-trigger efficiency for minimum-bias (solid 
crosses), PYTHIA charm (dashed crosses), and vrvr events (stars) vs. 

threshold in photoelectrons. 



estimated in Table |^ assuming (for the sake of illustration) that the optical 
trigger is used for all-hadronic modes (but not for leptonic modes, for which 
the ffist-level trigger rate should be sufficiently low to be recorded directly) 
and performs as estimated above. Although due to off-line selection cuts not 
yet simulated, realistic yields could be a factor ^2 — 3 below those indicated, 
the total reconstructed sample is well in excess of 10^ events. Given the factor 
^ 2 mass- resolution improvement compared to E791, one can infer a factor 
~ 50 improvement in statistical significance in a typical decay mode. Since 
the charm cross section at 120 GeV proton-beam energy may be several % of 
that at 800 GeV, and the geometrical acceptance remains ~ 50%, interesting 
charm sensitivity may also be available using Main Injector beam during 
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Table 2: Estimated yields of reconstructed events (antiparticles included) 



mode 


charm frac. 


BR 


acceptance 


efficiency 


yield 


D'^ Ktt 


0.6 


0.0365 


0.56 


0.1 


1.2 X 10^ 


D+ K*ny 


0.3 


0.027 


0.4 


0.25 


8 X 10^ 


— > KTr/ii/ 












all 


1 


0.1 


fa 0.4 


0.1 


4 X 10^ 



Tevatron Collider running; at the least, there will be opportunity to debug 
and test the spectrometer thoroughly so that full-energy beam may be used 
with optimal efficiency. 

4 Summary 

A fixed-target hadroproduction experiment capable of reconstructing in ex- 
cess of 10^ charm events is feasible using detector, trigger, and data acqui- 
sition technologies which exist or are under development. A typical factor 
~ 50 in statistical significance of signals may be expected compared to E687 
and E791. We expect the spectrometer sketched here to cost substantially 
less than HERA-S (whose cost was estimated at 33M DM in 1994). We 
anticipate an exciting future for charm physics at the turn of the century. 
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